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Nanoscale systems of metal atoms antiferromagnetically exchange coupled to several magnetic 
impurities are shown to exhibit an unconventional re-entrant competition between Kondo screening 
and indirect magnetic exchange interaction. Depending on the atomic positions of the magnetic mo- 
ments, the total ground-state spin deviates from predictions of standard Ruderman-Kittel-Kasuya- 
Yosida perturbation theory. The effect shows up on an energy scale larger than the level width 
induced by the coupling to the environment and is experimentally verifiable by studying magnetic 
field dependencies. 
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Introduction. The competition between different 
mechanisms working at comparable energy scales is in 
many cases the origin of unconventional physical prop- 
erties. The Kondo temperature Tk is the energy scale 
below which the local magnetic moment of a magnetic 
impurity is screened by the conduction-electron spins 
of a metallic host PQ. In multi-impurity systems, the 
Kondo effect competes with inter-impurity magnetic cor- 
relations caused by the Ruderman-Kittel-Kasuya-Yosida 
indirect magnetic exchange interaction [5H1]. As was al- 
ready pointed out by Doniach in 1977 [Jj, their different 
scaling with the antiferromagnetic coupling J, namely 
Tk ~ e _1 / J and Jrkky ~ J 2 (for weak J), gives rise to 
a point J — Jd where Tk = Jrkky- This roughly marks 
a crossover or even a phase transition and represents an 
important key to understand the phase diagrams of di- 
lute Kondo systems or Kondo lattices. Here we ask: How 
does this competition between the Kondo screening and 
RKKY interaction change if the quantum system is made 
so small that its conduction-electron spectrum becomes 
discrete? How does the physics change due to the pres- 
ence of a third energy scale, the level spacing A close to 
the Fermi energy? 

For a single magnetic impurity, i.e. for the "Kondo- 
box" problem [6] , there is already a competition, namely 
between A and T K [THE] : If the level spacing A becomes 
comparable to the bulk Tk, logarithmic Kondo correla- 
tions are cut, and the extension of the Kondo screening 
cloud is actually given by the system size. There is even- 
tually only a single conduction-electron state within the 
Kondo scale Tk around the Fermi energy which is avail- 
able to form the "Kondo" singlet. Here, wc will argue 
that this feature results in an unconventional, spatially 
dependent competition with the RKKY interaction for 
the multi-impurity case. This becomes relevant for mag- 
netic nanostructure physics in an important parameter 
range and for the bottom-up construction of spintronics 
devices. 

Experimentally, Kondo boxes have been realized as an 
individual grain [16l [17] , as a single metallocene molecule 
[18] , or as a small quantum dot which acts as a spin-half 
impurity and which is coupled to a large quantum dot 



with a finite level spacing [T5]. The Tk/A ratio can be 
tuned by varying the voltage on the gates separating the 
two dots. One-dimensional Kondo boxes can be realized 
by small Co clusters in short carbon nanotubes pieces [2U] 
or by a carbon nanotubc exchange coupled to a localized 
electron |21) . 

With recent progress in manipulation and character- 
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FIG. 1: (Color online) Competition between Kondo screen- 
ing and RKKY interaction in a quantum box. Local Kondo 
singlets (LKS) are formed for extremely strong J (regime 4). 
With decreasing J, the energy to break up a Kondo singlet 
becomes exponentially small (3). Below Jd, RKKY coupling 
is dominant (2). Finite-size effects set in for J < J a defined 
as the point where the bulk Kondo temperature equals the 
finite-size gap, Tk = A. If the system is sufficiently large, 
we have J a < Jd- The singlet-formation energy is linear in 
J for "good" and vanishes for "bad" sites (1): I.e. for sites 
where the kp conduction-electron wave function has a finite 
(vanishing) weight, the Kondo effect "wins" (RKKY wins). 
The remaining unscreened moments are subjected to nonlocal 
RKKY exchange subsequently. At a still lower energy scale 
S, determined by the residual coupling of the box to the envi- 
ronment, the Kondo scale becomes exponentially small again 
with a screening cloud leaking out into the environment (0) 
[9]. Here, for J < Js, the RKKY coupling is dominant again 
[141 115j . The arrow marks the range covered by our numerical 
calculations. 
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ization of magnetic systems on mesoscopic and nano 
scales, studies of the competition between Kondo screen- 
ing and RKKY interaction in a quantum box come into 
reach as, e.g., in artificial and tunable double quantum- 
dot systems [22] . Scanning-tunneling techniques nowa- 
days allow to measure atomically precise maps of the 
RKKY coupling between individual adatom pairs on 
metallic surfaces [23j and the bottom-up construction 
of tailored magnetic nanostructures with atomic control 

nam]. 

Main results. Here we show that there is an uncon- 
ventional competition between Kondo and RKKY that 
sets in for J < Jd in a system with several magnetic 
impurities coupled locally via J to the sites of a finite 
quantum box (see Fig. [lj. This re-entrant competition is 
a consequence of a combination of two important finite- 
size effects becoming vital if T]^ ulk < A or, equivalently, if 
J < J a- First, there is an odd-even effect in the number 
of conduction electrons N. For even N, the Fermi en- 
ergy lies in a finite-size gap ("off resonance"). As the 
Kondo scale is cut by A, the local moments are un- 
screened but couple via the RKKY interaction. For odd 
N ("on resonance"), however, cutting the Kondo scale 
implies that standard perturbation-theory in J is regu- 
larized and that already the first-order-in- J term leads 
to a singlet formation. This linear-va-J Kondo scale wins 
over the RKKY scale oc J 2 for weak J if, secondly, the 
single-conduction-electron eigenstate at the Fermi edge 
has finite weight at the respective impurity site ("good 
site" ) while for "bad sites" unscreened moments are sub- 
jected to RKKY coupling. We demonstrate that a fairly 
complete qualitative picture of this physics is obtained 
by studying t/iree-impurity Kondo models using exact 
diagonalization, perturbation theory and density-matrix 
renormalization group (DMRG) [26) I27j. 

For a completely isolated box the conduction-electron 
density of states consists of a set of delta-peaks separated 
by A in the vicinity of the Fermi edge. A weak coupling of 
the box to an environment (e.g. leads), however, turns the 
delta-peaks into resonances with a characteristic width 5 
[9j. Thereby a fourth scale comes into play. At this 
energy scale ~ S <C A, or equivalently for J ~ Jg <C J a, 
the novel re-entrant competition in the Kondo- vs. -RKKY 
quantum box breaks down. Our study thus bridges the 
gap between the conventional Kondo-vs.-RKKY physics 
for J > J a and the regime J < Js studied previously 

on eng. 

Multi-impurity model. We consider M spins S m , with 
spin-quantum numbers 1/2, which are coupled locally via 
an antiferromagnetic exchange J > to the local spins Sj 
of a system of N itinerant and non-interacting conduction 
electrons. An even total number N + M of conduction 
electrons and localized spins 1/2 is assumed such that a 
Fermi-liquid state with a total spin singlet can be reached 
for L —> oc. 

The conduction electrons hop with amplitude t = 1 



between non-degenerate orbitals and nearest-neighboring 
sites of a lattice with finite number of sites L. The Hamil- 
tonian is 

M 

n = -tJ2 C *W + J E S ^Sm ■ (1) 

{i,j),u m=l 

Here, Ci a annihilates an electron at site i — 1, ...,L with 
spin projection a =f, I, and Sj = 5 X^o-o-' c I<r CT c * ff/ is 
the local conduction-electron spin at i, where <x is the 
vector of Pauli matrices. Impurity spins couple to the 
local conduction-electron spins at the sites i m where m = 
1,...,M. 

For the sake of clarity, we consider a one-dimensional 
chain with open boundaries. Diagonalization of the tight- 
binding part in Eq. ([!]) yields nondegenerate conduction- 
band energies Sk = —2t cos k with discrete fc = nn/(L+l) 
labeled by integers n = 1, ...,L. For a half-filled system, 
i.e. for N — L electrons, this results in an energy-level 
spacing A = 2tsin(7r/(L+l)) at the Fermi energy ep = 
if L is odd. The local density of states at i consists of 
a finite number of <5-peaks only, pu(oj) = J2k ^ 2 /c^( w ~ 
(efc — £p)) where Uik — y/2/L sin(i k). Opposed to the 
continuum limit L — > 00, it is no longer finite in the 
vicinity of £f- 

Effective RKKY model. For two impurity spins and 
even TV, £p lies in a finite-size gap of the single- 
conduction-electron spectrum ( "off- resonant case"). 
There is a crossover, for finite L, from local Kondo- 
singlet formation at strong J to nonlocal RKKY cou- 
pling of the spins for J — > as can be seen in spin- 
correlation functions and susceptibilities [2"5l |2"5] . Note 
that for J —> (for Tk < A) the Kondo effect is absent, 
and free moments are generated: The impurity spins can- 
not dynamically couple to the conduction-electron sys- 
tem as the Fermi sea is non-degenerate and thus a finite 
energy ~ A would be necessary to screen a spin. Hence, 
the low-energy sector is exactly described by the effec- 
tive RKKY two-spin model -Hrkky = — J12S1S2 with 
J12 oc (~l)l 4 i-^l p/\i x - i 2 \ for J -> 0. This standard 
RKKY physics is easily recovered for very small L by full 
diagonalization of H. 

M = 3 spins. The quantum box with three impuri- 
ties is qualitatively different since a non-trivial competi- 
tion between RKKY coupling and Kondo-singlet forma- 
tion is possible even for J — >• 0: As N is odd, the high- 
est one-particle eigenenergy Sk F is singly occupied, and 
thus the ground state of the conduction-electron system 
is two- fold Kramers degenerate ("resonant case"). Hence, 
screening of an impurity spin is possible for J — > but 
competes with the RKKY exchange: Consider a chain 
with one impurity spin coupled to the system via J at 
the central site while the other two spins are coupled 
to the adjacent sites. At second order in J, perturba- 
tion theory predicts indirect antiferromagnetic coupling 
of the central with the adjacent spins and ferromagnetic 
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coupling between the latter. A naive argument based 
on Doniach's idea would first consider the correspond- 
ing effective three-spin RKKY model which has a to- 
tal impurity-spin doublet ground state (for any absolute 
magnitudes of the RKKY couplings). 

However, it is in fact the Kondo effect that wins in 
this case: Perturbation theory for the Kondo problem 
is regularized due to the finite-size gap A > and pre- 
dicts that, if an impurity spin dynamically couples to the 
conduction- electron system this happens on a linear-in- J 
scale. For sufficiently weak but finite J, this is larger than 
the RKKY scale oc J 2 and we thus expect formation of a 
spin singlet involving the conduction-electron spins. We 
will call this a "Kondo singlet" although 7k oc J rather 
than being exponentially small. 

Whether or not there is a perturbative coupling of the 
impurity spin, depends on the weight factor Uik at k = 
hp: Uik F is just the z-component of the conduction-band 
one-particle energy eigenstate at Ef- At hp = tt/2, we 
have Uik F = \/2/L sin(i kp ) 7^ for i — 1 , 3, L — 2, L. 
We call these sites "good", while Uik F = for "bad" 
sites i = 2,4, ...,L — 1. Assume that L = An + 3 with 
integer n, i.e. the central site is bad. The physics of 
this system is thus dominated by a Kondo screening of 
one of the two outer impurities. The RKKY interac- 
tion comes into play in a second step only and mediates 
an antiferromagnetic coupling between the two remain- 
ing spins which form a spin singlet. Hence, the ground 
state is a Kondo singlet entangled with an RKKY sin- 
glet, |GS) = |Ki> <g> |RKKY 23 ) - |RKKY 12 ) ® |K 3 ), and 
the total spin S to t = 0. 

Low- energy model. This has in fact been verified by 
exact diagonalization for very small systems with L = 3 
and L = 7 but also for somewhat larger chains with L = 
19 and L — 51 using DMRG (see Ref. [29] for methodical 
details). Moreover, the ground-state symmetry and all 
spin-spin correlation functions are, up to order J 2 , fully 
determined by an effective four-spin model (for M = 3) 
that replaces the RKKY Hamiltonian: 

M M 

H c ff = ^2 (Jm) + Jm^)Srn,SF ~ ^ JmnS m S n . (2) 

m— 1 m,n=l 

Here, the effective coupling constants depend on Ui m k 
and £k- sp is the spin of the fully delocalized /c^-electron 
- the Kondo cloud extends over the entire system for 
J — > 0. The ground-state symmetry crucially depends 
on the position of the impurities as both the linear, 
Jm — J\Ui m k F \ 2 , and the quadratic Kondo coupling, 

Jm = J 2 \Ui m k F \ 2 J2 P> k F ej-ek ' vanisn at sites where 
\U lmkF \ 2 = while the RKKY coupling J mn oc J 2 ^ 0. 

In a situation where the three local spins are all in- 
teracting with bad sites, for example, there is no Kondo 
effect. The inter- impurity distance is even which implies 
a ferromagnetic RKKY coupling leading to a total im- 
purity spin =3/2 that does not couple to the total 
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FIG. 2: (Color online) Inter-impurity and local spin correla- 
tions showing re-entrant competition between Kondo screen- 
ing and RKKY coupling for the M = 3 model ("bad-good- 
bad") with increasing (see arrows) system size up to L — 201 
sites. Note the log-scale for J. With decreasing J, different 
regimes 1-4 (cf. Fig. [T]) are found: (4) Strong-coupling limit 
with local Kondo singlets. (3) Kondo screening dominates 
over RKKY interaction. (2) RKKY interaction is the lead- 
ing energy scale for large systems (Jd > J a)- (1) The linear 
Kondo scale dominates and leads to a screening of the central 
spin (coupled to a good site, green) while in a second step 
RKKY couples the remaining spins (at bad sites, red) to a 
nonlocal triplet. 

conduction-electron spin Scond = 1/2 for weak J. Hence, 
the total spin S^ot = 1 or S tot = 2, the ground state is 
8-fold degenerate. 

If, at the same inter-impurity distance, all spins inter- 
act with good sites, the Kondo scale is large as compared 
to RKKY. Hence, one of the spins is Kondo-screened 
while the remaining ones are RKKY-coupled to a triplet: 
Stot = 1- The ground state is an entangled Kondo- 
singlet-RKKY-triplet state. 

Re-entrant competition. What happens for finite A 
as J increases and what happens for finite J if A de- 
creases, i.e. if the system size L increases? Let us discuss 
those questions for the remaining case "bad-good-bad" : 
Two spins Si and S3 couple to bad sites neighboring 
the central site of a chain with L = 4n + 1 . The central 
site is good. Therefore, for small J and small L we have 
y(buik) ^ ^ perturbative arguments given above ap- 
ply, the Kondo scale is linear in J, and S2 is Kondo 
screened. The weaker (ferromagnetic) RKKY interac- 
tion then couples S\ and S3 to a nonlocal spin triplet, 
and thus 5 t ot = 1- 
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This is nicely reflected in the different ground-state 
spin-correlation functions obtained by DMRG which are 
shown in Fig. [2j For L = 9 at J ->■ 0, we find (SiS 3 ) -> 
1/4 while (S 1 S 2 ) -> (and (S 2 S 3 ) -> by symmetry). 
The local correlations vanish (Sis^) — > at the bad sites 
but for the good one (S^s^) remains finite as J — > 0. 

The strong-J limit is also easily understood: For J > 
Ja the distinction between good and bad sites becomes 
irrelevant, and as J — > oo all local spin correlations tend 
to —3/4. This indicates local Kondo-singlet formation 
which is basically independent of the system size. 

At intermediate J, i.e. Ja < J < Jd, the conven- 
tional competition between Kondo screening and RKKY 
interaction should be recovered (since Ja < J) , and the 
RKKY interaction should win (since J < Jd). However, 
this crucially depends on the system size: For small sys- 
tems, we rather have Ja > Jd, and the intermediate- J 
regime is skipped. This can be seen in (S1S2) which 
stays close to zero in the entire J regime with an only 
shallow minimum around J = 1. 

Only if the system is sufficiently large, namely if 
Ja < Jd (i-e. A < Jrkky), RKKY correlations among 
the impurity spins can develop (this is also the situa- 
tion sketched in Fig. [lj. In fact, as is seen in Fig. [2] 
for L = 201, the correlation (S1S2) -> -1/2 in the 
intermediate-J regime and (S1S3) close to 1/4. These 
are exactly the spin correlations of a three-spin system 
with ferromagnetic coupling between Si and S3 and an- 
tifcrromagnetic ones else. With decreasing J, the system 
eventually crosses over to the perturbative regime at a 
Ja that strongly decreases with decreasing L. 

One might speculate that the nonlocal impurity spin 
doublet that is formed by the RKKY interaction in the 
intermediate-J regime, A < T^ ulk) < Jrkky, is Kondo 
screened in a subsequent step on a very small energy 
scale y( 3spiIls ) <g; T^ ulk \ A conventional Kondo effect 
would be obtained only if A < T^ spins) < T^ ulk) < 
Jrkky j i-e. for very large systems. Here, we rather expect 
2^3 spins) ^ ^ < jip>ulk) < Jrkky > i-e. the corresponding 
Kondo correlations are again cut by the system size and 
a linear-in-J 3-spin Kondo scale develops. Still, for the 
maximum system size considered here (L = 201), this 
scale appears too small to be resolved numerically. 

Another issue for the quantum-box regime j^, bulk ) < 
A, is a weak coupling to the environment. Eventually, 
the total spin S tot = 1 of the "bad-good-bad" configu- 
ration, for example, may be screened by environmental 
spin degrees of freedom. If the level broadening is small, 
S -C A, this is expected to take place on a very small 
energy scale <§C 5 [HI fl5j . In the "on resonance" case, 
however, the RKKY interaction is particularly enhanced 
by intra-resonance processes [15j . 

Conclusions. For an isolated quantum box including 
several localized magnetic moments coupled antiferro- 
magnetically to conduction-electron spins, there is a sub- 



tle competition between Kondo screening and RKKY in- 
teraction even in the limit J — > 0. The symmetry of 
the ground state crucially depends on the geometrical 
position of the impurities. It is predictable by perturba- 
tive techniques and experimentally accessible by study- 
ing field dependencies [35] . For sufficiently large systems 
( Jd > Ja), an increasing J drives the system from Kondo 
screening on a linear-in-J scale to conventional RKKY- 
coupled moments until the Kondo effect takes over again 
for large J. The relevant parameter range can be esti- 
mated roughly by setting Tk ~ A ~ Jrkky : RKKY cou- 
plings in the range of 0.1-10 meV correspond to Kondo 
temperatures of 1-100 K and, in a free-electron model [6], 
to system volumes of (11.5nm) 3 — (2.5nm) 3 . 

As a possible application, we suggest to employ 
quantum-confined systems at surfaces by means of 
scanning-tunneling techniques with the objective to con- 
struct nano-spintronics devices bottom- up |24j . Even if 
the confinement, e.g. due to insulating spacers, is not per- 
fect, it can lead to strong variations in the local density 
of states as is known from quantum corrals, for exam- 
ple [3D]. The corresponding strong spatial dependence 
of the Kondo temperature is sufficient, if combined with 
an atomically precise positioning of atoms, to utilize the 
Kondo-vs.-RKKY physics in a quantum box. 

The latter is limited for strong J by local Kondo-singlct 
formation and by residual couplings to the environment 
for very weak J. Generalizations to more impurities, sys- 
tems off half-filling and multi-orbital systems are worth 
to be explored. We also expect our work to serve as a 
bottom-up route to address dense or Kondo-lattice mod- 
els. 
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